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Abstract 

Quantum fluctuations in the density of a fluid with a linear phonon dispersion relation are stud- 
ied. In particular, we treat the changes in these fluctuations due to non-classical states of phonons 



qh| and to the presence of boundaries. These effects are analogous to similar effects in relativistic 



quantum field theory, and we argue that the case of the fluid is a useful analog model for effects 
in field theory. We further argue that the changes in the mean squared density are in principle 
observable by light scattering experiments. 
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I. INTRODUCTION 



It is well known that quantized sound waves, whose excitations are phonons, share sev- 
eral properties with relativistic quantum fields, such as the electromagnetic field. This is 
especially true when the phonon dispersion relation is approximately linear, which will be 
assumed throughout this paper. There is a phononic analog of the usual Casimir effect, but 
it tends to be quite small. For example, the force on two parallel plates due to phonon zero 
point energy is smaller than that in the electromagnetic case by the ratio of the speed of 
sound in the fluid to the speed of light [1]. This ratio is typically of order 10^^. However, 
forces due to classical stochastic sound fluctuations have been discussed recently by several 
authors [1, 0, 0, 0, 0, 0] , and can be larger. Here we will study the local changes in density 
fluctuations of a fluid due either to the presence of boundaries or to changes in the quantum 
state of the phonons. This is an analog of the effect of boundaries on the quadratic expecta- 
tion values of relativistic quantum flelds, such as the mean squared electric fleld. In a related 
context, Unruh [sl has shown that the velocity potential (J), of a moving fluid with velocity 
V = V0, satisfles the same equation as does a relativistic scalar fleld in a curved spacetime. 
The present paper is an expanded version of Ref. [9!]. We begin in Sect. [TTlby reviewing the 
quantization of sound waves in a fluid and the calculation of the density correlation function. 
We also review recent work on the scattering of light by zero point fluctuations in a fluid. In 
Sect. Illlt we consider the effects of a squeezed state of phonons on the local density fluctu- 
ations. The effects of boundaries are treated in Sect. llVt where several different geometries 
are discussed. Our results are summarized and discussed in Sect. El 



II. PHONONS AND DENSITY FLUCTUATIONS IN A FLUID 



A. Quantization and the Density Correlation Function 

We consider the quantization of sound waves in a fluid with a linear dispersion relation, 
= cs q, where Qg is the phonon angular frequency, q is the magnitude of the wave 
vector, and cs is the speed of sound in the fluid. This should be a good approximation 
for wavelengths much longer than the interatomic separation. Let po be the mean mass 
density of the fluid. Then the variation in density around this mean value is represented by 
a quantum operator, p(x., t), which may be expanded in terms of phonon annihilation and 



creation operators as 10| 



p(x,t)=5:(Vq+&i/*)' (1) 

q 

where 



Here is a quantization volume. The normalization factor in Eq. ([2D can be flxed by 
requiring that the zero point energy of each mode be ^hflq and using the expression for the 
energy density in a sound wave, 

U=^p\ (3) 
Po 
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In the limit in which oo, we may write the density correlation function as 

npo_ 



(p(x,t)p(x',t')) = J rf^gfi.e^^'^-^''-^^^*) , (4) 



where Ax = x — x' and At = t — t'. The integral may be evaluated to write the coordinate 
space correlation function as 

(,(x.i)p(x'.0) = -^ (Ax^-34W - 

This is of the same form as the correlation function for the time derivative of a massless 
scalar field in relativistic quantum field theory, {(p{:x.,t) iphd,t')) . (This analogy has been 
noted previously in the hterature. See, for example, Ref. Apart from a factor of po, 

these two quantities may be obtained from one another by interchanging the speed of light 
c and the speed of sound cs- If c ^ cs, then 

(0(x, t) 0(x', t')) ^ po (p(x, t) p(x', t')) . (6) 

In the limit of equal times, the density correlation function becomes 

Thus the density fluctuations increase as |Ax| decreases. Of course, the continuum descrip- 
tion of the fluid and the linear dispersion relation both fail as | Ax| approaches the interatomic 
separation. Also note the minus sign in Eq. ([7j). This implies that density fluctuations at 
different locations at equal times are anticorrelated. By contrast, when C5|At| > |Ax|, then 
(p(x, t) p(x',t)) > and the fluctuations are positively correlated. This is complete analogy 
with the situation in the relativistic theory. Fluctuations inside the lightcone can propagate 
causally and tend to be positively correlated. Fluctuations in a fluid for which C5I At| < | Ax| 
cannot have propagated from one point to the other, and are anti-correlated. This can be 
understood physically because an over density of fluid at one point in space requires an 
under density at a nearby point. 



B. Light Scattering by Density Fluctuations 

In Ref. [13] , the cross section for the scattering of light by the zero point density fluctu- 
ations is computed for the case that the incident light angular frequency is large compared 
to the typical phonon frequency. The result is 

f^) = 7^(13^ ^^^(e,,.e„,,)^ (8) 

where 9 is the scattering angle, V is the scattering volume, and t] is the mean index of 
refraction of the fluid. In addition, ek,A and ek',A' are the initial and final polarization 
vectors, respectively. The uj^ dependence of the scattering cross section can be viewed as 
the product of the cu^ dependence of Rayleigh-Brillouin scattering and one power of uj coming 
from the spectrum of zero point fluctuations in the fluid. The factor of r}^ represents the 
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influence of tlie fluid on light propagation before and after the scattering process, and arises 
as a product of a factor of rj in the incident flux and a factor of t]^ in the density of final 
states Because light travels through the fluid at speeds much greater than the sound 
speed, light scattering reveals a nearly static distribution of density fluctuations. Thus we 
can regard Eq. ([H]) as a probe of the fluctuations described by Eq. ([7]). The scattering by zero 
point fluctuations is inelastic, with the creation of a phonon. Thus, the scattering described 
by Eq. ([H]) is strictly Brillouin rather than Rayleigh scattering. 

This scattering by zero point density fluctuations should be compared to the effects of 
thermal density fluctuations. The ratio of the zero point to the thermal scattering for the 
Stokes line may be expressed as 



{d(T/dQ)TB ^ \2kBTj \ c ' ' 




(9) 



The index of refraction, rj, and the quantity po {dG/dpo)g, which involves a derivative of the 
fluid dielectric function with respect to density at constant entropy are both of order unity. 
Hence R is primarily determined by the ratio of the photon energy to the thermal energy, 
and the ratio of the speed of sound to the speed of light. 

Note that the zero point scattering cross section, Eq. ([H]), is the sole cross section at 
zero temperature. At finite temperature. Stokes line cross section (describing the process in 
which a phonon is emitted) is modified by the factor 

where (uq) is the mean number of phonons in mode q, and hflg is the phonon energy. In 
the low temperature limit, kT <C hQq, this correction factor goes to unity, giving the zero 
point result. In the high temperature limit, kT ^ hQg, it becomes 

k-T 1 

The leading term is the usual high temperature limit. The next term is the zero point effect, 
giving rise to a contribution to the cross section proportional to uj^. More precisely, it is 
1/2 of the zero point effect, the other 1/2 having been canceled by the thermal correction. 
Our view is that zero point fluctuations are always present at all temperatures, but in this 
case the thermal correction partially masks the zero point effect. However, the half which 
remains is potentially observable. For experiments in the high temperature limit, the lj^ 
part of the cross section is given by 1/2 of the right hand side of Eq. ([8]). 

Some numerical estimates for various fluids are given in Ref. [ij] for violet light with a 
wavelength of A = 350nm. For the case of liquid neon, R ^ 0.13, so that about 13% of the 
Stokes line is due to zero point motion effects [l3[ , which might be detectable experimentally. 
Even in the case of water at room temperature, R ~ 0.004. Although small in absolute terms, 
this is surprizingly large for a macroscopic quantum effect at room temperature. 

It is interesting to note that if one were to look only at the total Brillouin cross section 
(Stokes plus anti-Stokes), the zero point effect would be masked at high temperatures. The 
anti-Stokes line describes phonon absorption, so in the limit that u ^ Qq, its cross section 
is of the same form as that for the Stokes line, but its thermal correction factor is (nq). The 



4 



total cross section from both lines has a factor of 



2 (V + 1 = coth f ~ ^ + 0(1/T) , kT:$>hn,. (12) 



2kT J hVtq 



Here the thermal part completely masks the zero point part, leaving a residue of order 1/T. 
The same masking effect also occurs for the energy of a collection of harmonic oscillators, 
which is proportional to the quantity in Eq. (JT2l) . The thermal effect on scattering is often 



described by a structure factor. See, for example, Ref [IJ]. The hyperbolic cotangent form 
of the srtucture factor, corresponding to Eq. p^ . was calculated in Ref. (l5| . 

In the remainder of this paper, we will discuss modifications to the local density fluctua- 
tions due to the phonon state or to boundaries. These modifications are at least in principle 
observable through changes in the scattering cross section, Eq. ([8]). 



III. SQUEEZED STATES OF PHONONS 

Here we consider the case where the phonon field is not in the vacuum state, but rather a 
squeezed state. The squeezed states are a two complex parameter family of states in which 
the quantum uncertainty in one variable can be reduced with a corresponding increase in 
the uncertainty of the conjugate variable. See, for example, Refs. |16l . Il7l | for a detailed 
treatment of the properties of the squeezed states. We will focus attention on the case of the 
squeezed vacuum states \C) for a single mode, labeled by a single complex squeeze parameter 

C = re^^ (13) 

and defined by 

|C) = ^(C)|0). (14) 

Here 

5(C) = e^K*'^'-^("^)'l (15) 

is the squeeze operator and a and are phonon annihilation and creation operators for the 
selected mode. This set of states is of special interest because they are the states generated 
by quantum particle creation processes, and they can exhibit local negative energy densities. 



(See, for example, Refs. |l8l. Il9l]) 



Consider the shift in the mean squared density fluctuations between the given state and 
the vacuum 

(p')/?=(C|p'IC)-(0|p^|0), (16) 

the "renormalized" mean squared density fluctuation. The result for this quantity in a single 
mode squeezed vacuum state for a plane wave in the 2;-direction is 



(p )r = — ^— sinhr {sinhr — coshr cos[2{kz — ut) + S]} . (17) 



Here we have used the identities 

S\C)aS{0 = a cosh r-a^e^-^ sinhr, (18) 

and 

S^C)a^ S{C) = coshr -ae"*^ sinhr. (19) 
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Note that this quantity can be either positive or negative, but its time or space average is 
positive. The suppression of the local density fluctuations in a squeezed state is analogous to 
the creation of negative energy densities for a massless, relativistic field. (Compare Eq. (fT7|) 
with Eq. (48) and Fig. 8 in Ref. [T^.) 



IV. BOUNDARIES 

If we introduce an impenetrable boundary into the fluid, the phonon field will satisfy 
Neumann boundary conditions 

n ■ V5p = (20) 

as a consequence of the impenetrability. Thus there will be a Casimir force on the boundaries 
which is analogous to the Casimir force produced by electromagnetic vacuum effects. For 
example, consider two parallel plates, which will experience an attractive force per unit area 
of 

- = (21) 
A 480a4 ' ^ ' 

which is smaller than the electromagnetic case for perfect plates by a factor of cs/ (2c), and 
is thus quite small in any realistic situation. 

Henceforth, we consider the local effect of boundaries on mean squared density fluctu- 
ations, and now define {p^)r to be the change due to the presence of the boundary. This 
quantity is of interest both as an analog model for the effects of boundaries in quantum 
field theory, and in its own right. The shifts in density fluctuations are at least in principle 
observable in light scattering experiments. 

Our interest in the phononic analog model is inspired by the fact that the study of 
boundary effects in quantum field theory is an active area of research, and has given rise 



to some recent controversies in the literature [20|, |21|. One question is the nature of the 
physical cutoff which prevents singularities at the boundary. An example of the subtleties is 
afforded by the mean squared electric and magnetic fields near a dielectric interface. When 
the material is a perfect conductor, these quantities are proportional to z~'^, where z is the 
distance to the interface. Specifically, in Lorentz-Heaviside units their asymptotic forms are 

She 1 , , 

^ ) o-^ 22 

and 

B')--^A- (23) 

One might expect that a realistic frequency dependent dielectric function would remove this 
singularity, but this is not the case. Instead one finds [23] that 

E^) (24) 
/ 327r z^ ^ ' 

and 

(25) 

^ / 967rc ^2 ' ^ ^ 
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where ujp is the plasma frequency of the material. Thus some physical cutoff other than 
dispersion is required. For realistic materials, it is likely to be surface roughness, but fluc- 
tuations in the position of the boundary can also serve as a cutoff [23;]. In a fluid, there is 
always a physical cutoff at the interatomic separation. 

In the remainder of this paper, we will analyze (p^)ij in different geometries. 



A. One or Two Parallel Plane Boundaries 

In both of these case, the renormalized density two-point function may be constructed 
by the method of images. First consider the case of a single plate located at z = 0. Let 
Go denote the density correlation function in the absence of a boundary. The two-point 
function which satisfies the boundary condition Eq. fl20|) on this boundary is 



G = Go{At, Axr, z - z') + Go{At, Axy, z + z') (26) 
where is in the direction transverse to the plate. The renormalized two-point function is 

GH = Go(At,AxT,^ + /). (27) 
The resulting shift in the mean squared density is 



R 



327r2 Cs 



< 



(28) 



where z is the distance to the boundary. For the case of two parallel planes, the correlation 
function is given by an infinite image sum: 



G= [Go{At, Axr, z- z' - 2an) + Go{At, Axt, z + z' - 2an)] , 



(29) 



where a is the plate separation. If we use the identity 



E 



we can obtain the result 

{p')r 



6 dx"^ 



E 



(n — xy 



n^d^ 
6 dx'^ 



csc^(7ra;) 



(30) 



96 Cs 



1 3 — 2 sin^(7r2;/a) 
15 sin^(7rz/a) 



(31) 



where z is the distance to one boundary. Note that (p^)_R for both of these cases is negative 
everywhere. In the absence of a physical cutoff, both of these expressions diverge as z~^ 
near the boundaries, just as do the squared electric and magnetic fields near a perfectly 
reflecting plane. In contrast to the force between two plates, Eq. (12T!) . the shift in mean 
squared density is inversely proportional to the speed of sound, {p^)r oc l/c^. This is a 
general feature of all shifts due to boundaries. 
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B. A Three-Dimensional Torus 



Here we consider a rectangular box with periodic boundary conditions in all three spatial 
directions, with periodicity lengths Li, L2 and L3. Thus the three-dimensional space has 
the topology of S*^ x 5*^ x 5*^. This is closely related to the geometry of a waveguide, where 



the fluctuations of a relativistic scalar field were discussed by Rodrigues and Svaiter [24 
As in the parallel plane case, an image sum method may be employed to write 

00 

G= J2 GoiAt,x-x' + iLi,y-y' + mL2,z- z' + nL3). (32) 

£,m,n=— 00 



This leads to the result 



^^'^"^ 2Jls £ (i'Ll + m^Ll + n^LlY ' ^^^^ 

Here the prime on the summation indices denotes that the i = m = n = term is omitted. 
In this case, (p^)/?, is a negative constant. 



C. A Wedge 



Consider two intersecting plane which are at an angle of a with respect to each other. 
Now consider a point inside of this wedge which is located at polar coordinates (r, 6), where 
r is the distance to the intersection line and 6 < a. This geometry was treated for the 
relativistic case by Candelas and Deutsch [25I], whose Eq. (5.39) yields 



g2 



X, X 



(0') = lim —G(^,^, 

' x'^x dt^ ^ ' 3r2 B'^B 

Here Gr[B, 6') = G{d, 6') - Go{e, 6'), where 

Go{6, 6' 



lim 1 + 



Gn(e,e') 



'9-9'' 



is the empty space two-point function, and 



CSC 



Gi9,9') 



CSC 



7r{9-9') 
2^ 



+ CSC 



7i(9 + 9') 



2a 



is the two-point function in the presence of the wedge. 

We may combine these results to find for the phononic case 



hpo 



(34) 



(35) 



(36) 



X 



{(tt - a)(7r + a) sin^ {7i9 / a)[{7r^ + iW) sm^{7r9/a) - SOvr^] + 457r^} . (37) 



Again, this quantity is negative everywhere. 
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D. A Cosmic String 



As is well known, the space surrounding a cosmic string is a conical space with a deficit 
angle a < 27r. Quantum field theory in this conical space has been discussed by many 
authors, beginning with Helliwell and Konkowski j26|, and is similar to the wedge problem 
discussed above. Eqs. f lM|) and fl5S]) hold for the cosmic string as well as the wedge. However, 
Eq. fl5B]) is replaced by 



CSC 



^'1 



a 



(38) 



which is equivalent to Eqs. (15) and (16) in Ref. 2^. At a distance r from the apex, we find 



R 



14407r2 cs 

which is also negative everywhere provided that a < 27r. 



(27r - a)(27r + a)(lla^ + 47r^ 



(39) 



E. Near the Focus of a Parabolic Mirror 

The quantization of the electromagnetic field in the presence of a parabolic mirror was 



discussed by us in Refs. [27|, l28| , where a geometric optics approximation was employed to 
find the mean squared fields near the focus. This treatment lead to the result that these 
quantities are singular at the focus, diverging as an inverse power of the distance a to the 
focus. This result holds both for parabolic cylinders and for parabolas of revolution, and 
basically arises from the interference term of multiply refiected rays with nearly the same 
optical path length. The geometry is illustrated in Fig. [H An incoming ray at an angle of 
Q refiects at an angle of Q' to reach the point P, which is a distance a from the focus F, as 
illustrated. The distance from the focus to the mirror itself is 6/2 ^ a. 
The relation between Q and Q' is given by 

Q = \fm. (40) 

where 

, sin^6'' sin(6'' - 7) , . , 

/(^ ) = {\-co^Q') = + ^ ) ^^^(^ " ^) • ('^1) 

Note that Q is defined somewhat differently than in Refs. 27, 2^, so that now has 
the opposite sign. There will be multiply reflected rays whenever different values of Q' are 
associated with the same value of /. The function /(6'') is plotted in Fig.[2]for various values 
of 7. We can see from these plots that in general there can be up to four reflected angles Q' 
for a given incident angle Q. However, if the mirror size is restricted to be less than ^n/?), 
then there will never be more than two values of 9' for a given 9. Throughout this paper, 
we will assume ^0 < 27r/3, and hence have at most two reflected rays for a given incident 
ray. The two reflected rays will occur at ^' = a and 9' = j3, where 

/(«) = /(/?) • (42) 
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FIG. 1: The geometry of rays reflecting from a parabolic mirror is illustrated. An incoming ray at 
an angle of 9 reflects at an angle of 6' to reach the point P, which is a distance a from the focus 
F, and at an angle of 7. 





















1 ■ 


e' 


'~~----3.4^J__,,--^.e -0.8 











-1 ■ 







FIG. 2: The function f{0') is plotted for various values of 7. This function relates the angle of the 
incident ray, 6, to the angle of the reflected ray, 6', through the relation 9 = (a/b) f{6'). 

The difference in the optical paths of these two rays {j3 path minus a path) is denoted by 



A£. The detailed expression for this distance A£ used in Refs. [27|, l28| is not quite correct, 
as was pointed out to us by Vuletic The expression used in Refs. [23, 28], which we 

will denote by A£i, is the difference in distance traveled by the two rays after they cross a 
line of constant x, perpendicular to the axis of the mirror. This difference is 

A£i = a I cos 7(cos a — cos [3) + sin 7(sin a — sin /?) | . (43) 

However, the difference in optical path lengths is the difference in distance travel-led after 
crossing a line perpendicular to the incoming rays, as illustrated in Fig. [3l and is 

M = Ail - A^a ■ (44) 
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FIG. 3: Two rays reflecting from a parabolic mirror to the point P are illustrated. The first arrives 
at an angle of 9' = a, and the second at 9' = [5. The points of intersection with the mirror are 
(xi, yi) and (x2, ^2), respectively. The difference in path lengths (lower path minus upper path) is 
= Ml - M2, where M2 is illustrated. 

The correction term, A£2, is 

/S.^2 = o, [sin P sin(/5 — 7) — sin a sin(a — 7)] . (45) 
The corrected expression for A£ is then 

Ai = a cos 7(cos a — cos P + sin^ a — sin^ /?) + sin 7(sin a — sin /5 + sin /? cos /? — sin a cos a) 

(46)" 

The mean squared electric field near the focus of a parabola of revolution is, in the geometric 
optic approximation, 

(^V^l^/^, (47) 
The corresponding expression for a parabolic cylinder is 

{E\c = ^{E\r. (48) 
Note that in Eq. fHTI) . the integration is over 6, the angle of the incident ray, not 6', the 



reflected angle, as was incorrectly stated in Refs. [27|, 128 . 

A detailed discussion of the electromagnetic case will be given elsewhere. Here we are 
concerned with {p'^)r for the parabola of revolution, which is obtained from Eq. fHTl) by 
letting c — > cs and dividing by 2, leading to the result 

^'^-=2^sLwr' '''' 

where the factor of h has been restored. The corresponding expression of a parabolic cylinder 
is obtained by multiplying by 16/157r. Although the integrand in the above expression is 
singular at Ai = 0, it may be treated as a distribution and the integral is well defined. 
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Here we will treat only the case 7 = 7r/2, where the integrations may be done in closed 
form. In this case, (3 = —a, as may be seen from the fact that f{6') is now an even function: 



f{e') = {i + cose')cos{e'). (50) 

The minimum value of 6 in Eq. (H9|) is 6min = (O'/b) f{6o), where 6q is the angular size of the 
mirror. The maximum value in our case is 6max = 2a /b, corresponding to a = 0. We have 
that 

— = - / (a) = — - sina (2 cosa + 1) . (51) 
da b b 

This relation may be used to express as 

3/i J sin a (2 cos a + 1) 



or as 



3h fOo , 2cosa + l 3^ f^o , 2cosa + l 



{p)r Z2'k'^ csa^b jo sin^ q; (1 — cos a)^ QAn'^ cs a? b i -Oq | sin^ a| (1 — cos a)^ 

(53) 

This integral may be performed explicitly, with the result 

"Ml 

{p')r = -u 9i0o) , (54) 

409o7r"^ Cs a o 

where 

/ 1 + cos ^0 \ 30 cos^ ^0 - 120 cos^ + 160 cos^ - 40 cos^ - 94 cos - 224 

~ °^l^l-cos^o;^ 15(l + cos^o)(l-cos^o)' ■ 

(55) 

The function g{6Q) is negative everywhere, and is plotted in Fig. HI The singularity as 
represents a breakdown of the geometric optics approximation, as diffraction effects become 
more important for small ^o- For fixed 6q, the result is of the form 

(p% = -^<0, (56) 
csba-^ 

where C is a constant small compared to unity. 

This, and the analogous expressions for {E"^) and (-B^), which also are proportional to 
l/(6a^), are striking in that they can be large when the focus is far from the mirror itself, 
b ^ a. This result is controversial, and seems to be in conflict with a general result by 
Fewster and Pfenning 33|, which implies that quantities such as {E"^) or (/0^)_r should be 



proportional to the inverse fourth power of the distance to the mirror, which is to say oc 
in this case. On the other hand, there is a simple physical argument to the contrary, which 
we find compelling: the interference term between multiply reflected rays is slowly oscillating 
when A£ oc a is small, and should give a contribution proportional to an inverse power of a, 
as in Eq. (156|) . In any case, the study of the phononic case provides an additional theoretical, 
and potentially experimental, probe to better understand this issue. 
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FIG. 4: The function g{6o), defined in Eq. ([55]) . is plotted. 



V. SUMMARY AND DISCUSSION 

In this paper, we have treated the effects of squeezed phonon states and of boundaries 
on the local quantum density fluctuations of a fluid, assuming a linear phonon dispersion 
relation. The purpose of this investigation is two-fold. The modified density fluctuations are 
of interest in their own right and are in principle observable by light or neutron scattering. 
Secondly, the phononic system studied here is a potentially useful analog model for better 
understanding quantum fluctuations in relativistic quantum field theory with boundaries. 
After reviewing the density fluctuations in a boundaryless system in the phonon vacuum 
state, we treated the effects of a squeezed vacuum state of phonons. Here we found that 
such a state will have both local increases and local decreases in the mean squared density. 
However, the time or spatial averaged effect is an increase. This is in complete analogy 
to the case in relativistic quantum field theory, with the decrease in mean squared density 
corresponding to regions of negative energy density. 

We next turned our attention to the effects of perfectly reflecting boundaries and studied 
the cases of one and two parallel plates, a torus, a wedge, a cosmic string, and a parabolic 
mirror. In all of the cases examined, we found a decrease in mean squared density, {p'^)r < 
0. This amounts to a suppression of the usual zero point fluctuations, and is analogous 
to the suppression of vacuum fluctuations which can lead to negative energy density in 
quantum field theory. In general, {p^)r due to boundaries is inversely proportional to the 
speed of sound, cs- This is in contrast to total energies or forces, such as Eq. fl?Il) . which 
are proportional to cs, and to the mean squared electric or magnetic fields near a perfect 
reflector, Eqs. (l22l) and (!23l) . which are proportional to the speed of light. 

The case of the parabolic mirror is of particular interest. Here we were able to correct 
certain aspects of our previous treatment [i^, 2^ for electromagnetic fields. We find that near 



the focus, {p )r grows as the inverse cube of the distance ot the focus. For the phononic case. 
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this growth necessarily stops as the scale of interatomic spacing is reached. However, the 
analysis performed here for phonons also applies to the case of the quantized electromagnetic 
field, where one expects the same rate of growth in the mean squared electric and magnetic 
fields. 
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